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ABSTRACT 

One  of  the  main  thrusts  in  current  Army  Science  &  Technology  (S&T)  activities  is  the 
development  of  occupant- centric  vehicle  structures  that  make  the  operation  of  the  vehicle  both 
comfortable  and  safe  for  the  soldiers.  Furthermore,  a  lighter  weight  vehicle  structure  is  an 
enabling  factor  for  faster  transport,  higher  mobility,  greater  fuel  conservation,  higher  payload, 
and  a  reduced  ground  footprint  of  supporting  forces.  Therefore,  a  key  design  challenge  is  to 
develop  lightweight  occupant-centric  vehicle  structures  that  can  provide  high  levels  of  protection 
against  explosive  threats.  In  this  paper,  concepts  for  using  materials,  damping  and  other 
mechanisms  to  design  structures  with  unique  dynamic  characteristics  for  mitigating  blast  loads 
are  investigated.  The  Dynamic  Response  Index  (DRI)  metric  [1]  is  employed  as  an  occupant 
injury  measure  for  determining  the  effectiveness  of  the  each  blast  mitigation  configuration  that  is 
considered.  A  model  of  the  TARDEC  Generic  V-Hull  structure  comprises  the  vehicle  considered 
in  this  study.  The  material  properties  and  the  configuration  of  the  inner  bulkheads  that  connect 
the  outer  V-shaped  bottom  with  the  inner  floor  are  used  as  design  parameters  for  reducing  the 
DRI  at  a  typical  occupant  location.  In  this  particular  example,  it  is  demonstrated  that  the  weight 
of  the  structure  can  be  reduced  by  about  -12%  and  simultaneously,  the  DRI  can  be  reduced  by 
-24%).  This  is  achieved  by  creating  an  energy  absorbing/ decoupling  mechanism  between  the 
outer  hull,  the  inner  floor,  and  the  single  degree  of freedom  upper  torso  system. 

INTRODUCTION 

One  of  the  main  thrusts  in  current  Army  Science  &  Technology  (S&T)  activities  is  the 
development  of  occupant-centric  vehicle  structures  that  make  the  operation  of  the  vehicle  both 
comfortable  and  safe  for  the  soldiers.  Furthermore,  a  lighter  weight  vehicle  structure  is  an 
enabling  factor  for  faster  transport,  higher  mobility,  greater  fuel  conservation,  higher  payload, 
and  a  reduced  ground  footprint  of  supporting  forces.  Therefore,  a  key  design  challenge  is  to 
develop  lightweight  occupant-centric  vehicle  structures  that  can  provide  high  levels  of  protection 
against  explosive  threats.  Full  system,  end-to-end  [8,13-16]  as  well  as  Reduced  Order  [17-19] 
Modeling  and  Simulation  methodologies  are  extensively  used  for  the  development  of  blast¬ 
worthy  ground  vehicles  in  the  Army  Acquisition  process. 
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The  study  presented  in  this  paper  investigates  through  simulation  how  the  properties  of 
materials  used  in  the  construction  of  a  vehicle  structure  can  alter  the  dynamic  characteristics, 
offering  improved  isolation  from  the  blast  loads  at  a  lower  structural  weight.  A  Generic  V-Hull 
structure  developed  by  TARDEC,  aka  the  TARDEC  Generic  Hull  [13]  comprises  the  numerical 
model  for  investigating  these  concepts.  The  DRI,  which  is  a  standard  occupant  injury  metric  [1] 
for  underbody  simulations  and  testing,  is  used  in  this  study  as  a  measure  of  the  structural 
performance  with  respect  to  Survivability.  In  the  absence  of  an  actual  Anthropomorphic  Test 
Device  (ATD)  in  numerical  models  to  measure  lumbar  loads,  the  DRI  is  the  next  best  indicator 
of  lumbar  injury  performance  [10],  and  can  be  easily  calculated  from  structural  vehicle 
accelerations  as  shown  in  Appendix  A. 

In  the  literature,  various  concepts  of  employing  the  properties  of  materials  as  a  mechanism  to 
absorb  energy  have  been  presented.  For  example,  utilizing  shear  thickening  fluid  due  to  its  large 
capacity  for  energy  absorption  has  been  investigated  [2-3,  9].  Shear  thickening  fluid  is  a  specific 
type  of  non-Newtonian  fluid  with  its  viscosity  dependent  on  the  strain  rate.  It  acts  like  a  solid 
when  experiencing  a  large  shear  load,  such  as  an  impulse  of  high  pressure  but  of  short  duration 
from  a  blast,  and  returns  to  liquid  form  when  the  load  is  removed.  Employing  the  plastic 
deformation  induced  in  material  for  absorbing  energy  has  been  considered  for  designing  blast- 
resistant  structures  [4].  The  failure  mechanisms  in  unidirectional  fiber  reinforced  composites  of 
delamination,  fiber/matrix  debonding,  matrix  cracking,  and  fiber  breakage  have  been  considered 
for  creating  blast  mitigation  configurations  [5].  For  similar  purposes,  functional  graded  metallic 
materials  constructed  in  a  layered  sandwich  configuration  with  several  absorption  layers  have 
been  also  considered  [6]. 

The  concept  of  using  properties  of  “softer”  materials  is  investigated  in  this  paper.  It  allows  for 
higher  deformation  levels  in  the  structure  which  in  combination  with  a  damping  mechanism,  can 
result  in  a  reduced  base  excitation  leading  to  lower  DRIs  and  hence  occupant  injuries. 
Specifically,  the  properties  of  the  inner  bulkheads  that  connect  the  outer  V-Hull  bottom  to  the 
inner  floor  (Figure  3)  are  tuned  in  this  manner,  thereby  offering  an  isolation  mechanism  that 
reduces  the  DRI  metric.  In  the  following  sections  of  this  paper,  information  is  first  presented 
about  the  numerical  models  employed  in  this  study,  namely  the  TARDEC  V-Hull  model  and  the 
DRI  models.  The  software  code  LS-DYNA  is  used  in  the  blast  simulations  and  the  viscoelastic 
material  definition  is  used  for  defining  the  various  properties  of  the  internal  bulkheads  in  the 
parametric  study.  Therefore,  a  brief  discussion  on  the  viscoelastic  material  definition  in  LS- 
DYNA  is  presented.  Then,  two  different  lumped  parameter  models  for  the  DRI  metric  are 
described.  In  the  first  setup,  a  spring-mass  model  with  a  single  degree  of  freedom  (DOF) 
representing  the  upper  torso  of  the  occupant  (Appendix  A)  is  mounted  directly  in  the  middle  of 
the  inner  floor.  In  the  second  setup,  a  three-DOF  model  representing  the  upper  torso  of  the 
occupant,  an  energy-absorbing  seat,  and  an  energy-absorbing  floor  (Appendix  B)  is  mounted  to 
the  hull.  Finally,  the  process  followed  in  the  parametric  study  is  discussed  along  with  the  final 
design  configurations  that  reduce  both  the  structural  weight  and  the  DRI  metric. 

NUMERICAL  MODELS 

The  TARDEC  Generic  V-Hull  structure  is  presented  in  Figure  1.  It  comprises  a  representative 
but  generic/notional  vehicle  structure  that  can  be  used  in  survivability  research  studies.  The 
TARDEC  Generic  Hull  experiment  was  designed  to  provide  a  notional  geometry  for  underbody 
blast  analysis  and  to  evaluate  blast  mitigation  technologies.  Historically,  the  Department  of 
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Army  has  had  difficulty  collaborating  with  industry  and  academia  on  underbody  blast  events  due 
to  the  sensitive  nature  of  the  work.  Data  generated  from  testing  military  vehicles  is  typieally 
sensitive  and  not  readily  sharable.  To  alleviate  this  issue,  TARDEC  has  fabrieated  this  generic 
vehicle  hull  with  the  intent  to  share  data  with  academia  and  the  industry  to  spur  innovation  in 
blast  mitigation  teehnologies.The  main  dimensions  and  the  geometry  of  the  V-Hull  were  used  for 
ereating  a  simplified  model  that  is  used  in  this  work  (Figure  2).  The  simplified  model  has  the 
same  thieknesses  and  material  properties  for  the  main  struetural  eomponents  with  the  TARDEC 
V-Hull  strueture.  It  also  eontains  inner  bulkheads  conneeting  the  outer  V-shaped  bottom  with 
the  inner  floor  (Figure  3).  The  material  properties  of  the  bulkheads  are  used  as  design 
parameters  in  the  parametrie  study.  The  airblast  loading  feature  in  LS-DYNA 
(*LOAD_BLAST)  is  used  to  represent  a  mine  blast  threat  of  TNT  placed  0.2m  below  the  bottom 
eenter  of  the  vehiele. 


Figure  2.  Simplified  model  of  the  TARDEC  V-Hull  structure  that  is  used  in  this  study 

Figure  3a  depicts  the  single-DOF  (SDOF)  lumped  parameter  model  eonneeted  direetly  to  the 
vehiele  structure  to  evaluate  the  DRI.  The  upper  part  of  the  strueture  and  the  inner  floor  are 
removed  from  the  figure  in  order  for  the  internal  bulkheads  to  be  visible.  These  and  all  other 
parts  that  are  omitted  from  the  figure  for  visualization  purposes  are  included  in  the  simulations. 
The  DRI  is  used  as  a  metrie  for  assessing  the  safety  design  characteristics  of  a  vehicle.  It 
represents  the  dynamic  response  of  the  lower  lumbar  region  of  an  oceupant.  The  DRI  is 
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computed  from  the  maximum  dynamic  compression  measured  in  the  spring,  which  is  determined 
from  the  governing  equations  shown  in  Appendix  A. 

In  the  LS-DYNA  model,  a  single  spring-mass-damper  system  of  mass  34.51  kg,  with  natural 
frequency  equal  to  52.9  rad/s  and  damping  ratio  of  0.224  is  included  in  the  finite  element  model 
of  the  simplified  V-Hull.  As  mentioned  earlier,  a  second  three-DOF  lumped  parameter 
configuration  is  also  considered,  with  the  two  intermediate  DOFs  representing  the  seat  and  the 
floor.  This  second  configuration  is  presented  in  Figure  3b.  In  this  case,  the  DRI  is  determined 
by  the  relative  compression  in  the  spring  between  the  top  DOFs.  The  nonlinear  spring  constants 
for  the  lowest  (floor)  and  for  the  middle  (seat)  DOFs  are  shown  in  Figures  4a  and  4b, 
respectively,  and  a  damping  coefficient  of  9,220  (N*sec)/m  is  used  for  the  corresponding 
dampers.  By  including  these  one-DOF  and  three-DOF  systems  directly  in  the  simulations,  the 
calculations  in  the  governing  equations  shown  in  Appendix  A  and  B  are  automatically  performed 
by  LS-DYNA. 


Figure  3.  (a)  SDOF  model  (left);  (b)  Three-DOF  model  (right)  for  determination  of  DRI 


Nonlinear  Spring  Stiffness  for  Lower  sdofs  Nonlinear  Spring  Stiffness  for  Middie  sdofs 


Spring  Deformation  [m]  Spring  Deformation  [m] 


Figure  4.  Spring  stiffness  curves  for  the  floor  (left)  and  seat  (right)  DOFs 


The  density,  the  modulus  of  elasticity  and  the  dissipation  properties  of  the  material  comprising 
the  bulkheads  are  used  in  the  parametric  study  for  reducing  simultaneously  the  weight  of  the 
structure  and  the  DRI.  The  viscoelastic  material  definition  of  LS-DYNA  (MAT_061)  is  used  for 
modeling  this  material  [7],  which  models  both  viscous  and  elastic  characteristics  with  a  stress- 
strain  relation  that  depends  on  the  load  history.  It  behaves  as  a  spring-damper  system  and  two 
classical  models  (Maxwell’s  and  Kelvin’s)  are  employed  in  the  material  definition.  The 
parameters  which  are  considered  include:  mass  density,  bulk  modulus,  short-time  shear  modulus, 
long-time  shear  modulus,  and  a  decay  constant.  The  bulk  modulus,  the  short-time  shear 
modulus,  and  the  long-time  shear  modulus  are  determined  by  the  instantaneous  effective  spring 
coefficient  and  the  asymptotic  effective  spring  coefficient.  A  linear  relationship  between  the 
instantaneous  effective  spring  coefficient  Eq  and  the  asymptotic  effective  spring  coefficient  Eoo  is 
used,  Eo=1000*Eoo.  Therefore,  the  Poisson’s  ratio  of  the  material,  the  mass  density,  the  decay 
constant,  and  the  asymptotic  spring  coefficient  are  sufficient  for  defining  the  viscoelastic 
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material  properties.  In  this  study,  these  material  properties  are  used  for  creating  an  isolation 
mechanism  to  reduce  the  occupant  DRI  at  a  reduced  overall  weight. 

PARAMETRIC  STUDY 

The  configuration  with  the  single-DOF  model  cormected  to  the  floor  (Figure  3a)  for  evaluating 
the  DRI  was  analyzed  first.  Initially,  the  parametric  study  attempted  to  change  the  density  and 
the  stiffness  properties  of  the  entire  volume  of  each  bulkhead  part.  Two  main  conclusions  were 
drawn  from  this  initial  effort: 

•  First,  it  was  decided  to  preserve  the  original  steel  properties  for  the  upper  part  of  each 
bulkhead  (colored  yellow  in  Figure  5)  and  alter  the  stiffness,  the  density,  and  the  dissipation 
characteristics  in  the  remaining  portion  of  each  bulkhead  (colored  green  in  Figure  5).  The 
reason  for  this  approach,  is  to  avoid  excessive  local  flexibility  at  the  location  where  the 
SDOF  model  is  attached  to  the  floor  when  the  bulkhead  has  reduced  stiffness  properties.  The 
local  flexibility  at  the  attachment  point  makes  it  difficult  to  control  the  spring  compression 
that  determines  the  DRI. 

•  The  second  observation  was  that  the  overall  mass  of  the  vehicle  has  an  impact  to  the  overall 
rigid  body  response  of  the  vehicle  and  thereby  the  DRI.  Four  equal  lumped  masses  were 
added  at  the  four  comers  of  the  vehicle  to  keep  the  total  weight  constant  at  the  typical 
representative  weight  of  such  a  vehicle.  For  each  configuration,  the  values  of  the  lumped 
masses  were  selected  in  a  manner  that  the  overall  mass  of  the  vehicle  remained  constant. 

This  approach  also  reflects  the  ability  to  increase  the  payload  of  a  vehicle  even  when  the 
stracture  itself  becomes  lighter.  The  locations  where  the  lumped  masses  were  attached  are 
shown  in  Figure  5. 


Figure  5.  Partition  of  the  bulkheads  into  two  sections  (yellow  and  green);  and  locations  where 
lumped  masses  are  attached  for  preserving  the  overall  vehicle  mass 

The  final  configuration  identified  from  the  parametric  study  involving  numerous  configurations, 
had  the  following  values:  density  equal  to  6,000kg/m^,  asymptotic  spring  coefficient  equal  to 
3.8*10^  N/m,  and  a  decay  constant  equal  to  0.001 .  The  Poisson’s  ratio  did  not  vary  and  was  set 
equal  to  0.3.  Figure  6  summarizes  the  time  histories  of  deformation  of  the  upper  torso  relative  to 
the  pelvis,  for  the  original  configuration  (bulkheads  made  out  of  regular  steel)  and  the  final 
configuration.  The  values  for  the  maximum  spring  compression  and  the  associated  DRI  are  also 
included  in  the  Figure.  An  improvement  of  24.2%  is  observed  in  the  DRI  while  achieving  a 
12.5%  reduction  in  the  mass  of  the  structure  (as  explained  earlier  the  overall  mass  of  the  vehicle 
is  retained  constant  for  all  configurations).  Figures  7a  and  7b  present  the  displacement  histories 
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for  the  attachment  location  on  the  floor  where  the  SDOF  is  connected  (location  A),  and  for  the 
upper  torso  DOF  for  the  one-DOF  lumped  parameter  DRI  model  (location  B);  the  results  are 
shown  for  the  original  and  the  final  design  configurations,  respectively.  Being  total 
displacements,  these  results  include  the  effects  of  both  the  flexible  and  the  rigid  body  motion.  It 
can  be  observed  that  the  changes  in  the  material  properties  of  the  bulkhead  alter  primarily  the 
dynamic  response  of  the  connection  point  on  the  floor  and  in  this  manner  reduce  the  maximum 
relative  compression  and  the  DRI. 


0.015  ’ 


Outcome  from  parametric  study 
Regular  steel -original  configuration 


Figure  6.  Time  histories  of  deformation  of  the  upper  torso  relative  to  the  pelvis  in  the  one-DOF 

lumped  parameter  DRI  model 


Time 


Figure  7.  Dynamic  displacements  for  (a)  baseline  configuration  (top),  (b)  configuration  resulting 
from  parametric  study  (bottom).  Point  A  is  the  pelvis  (also  floor  attachment  location),  and  B  is 
the  upper  torso  DOF  of  the  one-DOF  lumped  parameter  model. 


In  a  similar  manner,  a  parametric  study  was  conducted,  involving  numerous  configurations,  for 
the  three-DOF  lumped  parameter  model  for  evaluating  the  DRI  (Figure  3b).  This  parametric 
study  retains  the  value  of  6,000  kg/m^  for  the  lower  section  of  each  bulkhead,  based  on  the 
configuration  identified  by  the  earlier  work.  The  asymptotic  stiffness  and  the  decay  constant 
comprise  the  varying  parameters.  For  the  final  configuration,  these  two  parameters  acquire 
values  of  2*10^  N/m  and  0.005,  respectively.  Figure  8  summarizes  the  time  histories  of 
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deformation  of  the  upper  torso  relative  to  the  pelvis,  for  the  original  configuration  (bulkheads 
made  out  of  regular  steel)  and  the  final  configuration.  The  values  for  the  maximum  spring 
compression  and  the  associated  DRI  are  also  included  in  the  figure.  The  improvement  in  the 
DRI  is  34.2%  this  time,  while  the  reduction  in  the  mass  of  the  structure  remains  at  12.5%. 
Figures  9a  and  9b  present  the  displacement  histories  for  the  attachment  location  on  the  hull 
where  the  three-DOF  lumped  mass  system  is  connected  (location  A),  and  for  the  upper  torso 
DOF  for  the  three-DOF  lumped  parameter  DRI  model  (location  B),  respectively.  In  this  case, 
both  the  response  of  the  upper  torso  DOF  and  the  connection  point  are  affected  by  the  material 
properties  of  the  bulkheads. 


Outcome  from  parametric  study 
Regular  steel -original  configuration 


Figure  8.  Time  histories  of  deformation  of  the  upper  torso  relative  to  the  pelvis  in  the  three-DOF 

lumped  parameter  DRI  model 


Figure  9.  Dynamic  displacements  for  (a)  baseline  configuration  (top),  (b)  configuration  resulting 
from  parametric  study  (bottom).  Point  A  is  the  pelvis,  and  B  is  the  upper  torso  DOF  of  the  three- 

DOF  lumped  parameter  model 
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CONCLUSIONS 

This  paper  demonstrates  how  material  properties  can  be  used  for  tuning  the  dynamic  behavior 
of  a  vehicle  to  simultaneously  reduce  the  weight  of  the  vehicle  structure  and  occupant  injuries. 
The  intent  is  not  to  identify  a  specific  material  or  design,  but  rather  exercise  a  process  for 
identifying  suitable  stiffness,  inertia,  and  damping/absorptive  characteristics  of  the  various 
components.  The  results  depend  on  how  and  where  the  seat  is  connected  to  the  vehicle,  the 
relative  stiffnesses  and  energy  absorption  characteristics  of  stroking  floors  and  seats  [1 1],  as  well 
as  the  material  properties.  The  selection  process  is  driven  by  controlling  and  minimizing  the 
energy  that  reaches  the  occupant  from  the  blast  and  the  resulting  occupant  injuries. 

ACRONYMS 


ARC 

Automotive  Research  Center 

ATD 

Anthropomorphic  Test  Device  (Dummy) 

CotS/COTS 

Commercial-Off-the-Shelf 

DA 

Department  of  the  Army 

DoB/DOB 

Depth  of  Burial 

DoD/DOD 

Department  of  Defense 

DOF 

Degree  of  Freedom 

DRI 

Dynamic  Response  Index 

DTIC 

Defense  Technical  Information  Center  http://www.dtic.mil 

EA 

Energy  Absorption 

FEA/FEM 

Finite  Element  Analysis/Model 

LS-DYNA 

COTS  structural  dynamics  software  from  Livermore  Software  Technology  Corporation,  CA 

M&S 

Modeling  &  Simulation 

MDOF 

Multiple  Degree  of  Freedom 

NA&ME 

Naval  Architecture  and  Mechanical  Engineering 

NATO 

North  Atlantic  Treaty  Organization 

PM 

Program  Manager 

R&D 

Research  &  Development 

RDECOM 

Research,  Development  and  Engineering  Command 

S&T 

Science  and  Technology 

SDOF 

Single  Degree  of  Freedom 

STANAG 

Standardization  Agreement 

TARDEC 

Tank  Automotive  Research,  Development  and  Engineering  Center 

UBB 

UnderBody  Blast 

UBM 

Underbody  Blast  Modeling/Methodology 
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Appendix  A:  Dynamic  Response  Index  (DRI)  -  SDOF  Mechanical  Model  [10] 


La 


•  6  is  the  relative  displacement 
between  the  upper  body  and 
pelvis  (6—A1-A2) 

•  ?  is  the  damping  coefficient^^ 
(0.224)  C  =  c/(2*sgn(m*k)) 

•  is  the  natural  frequency 
(52.9  rad/s)  %= 

•  Normalized  Lumbar  Force 

=  (k*5„  J  /  (mg)  = 

vWWVvV  s-Vs-Vs-Wv" 

dri  =  ■  " 

g 


m  =  34.51  kg 
k  =  9.66E04  N/m 
c  =  818.1  Nsec/m 

therefore: 

(D„  =  52.9  rad/s 
^  =  0.224 
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Appendix  B:  Dynamic  Response  Index  (DRI)  -  Three-DOF  Mechanical  Model  [11, 12] 


*  Takes  hull  acceleration  as  input. 

•  Accounts  for  energy  absorption  by 
the  floor  and  seat. 

•  Springs  representing  the  floor  and 
seat  are  piecewise-linear. 

*  The  spring  representing  the  spine  is 
still  linear. 


3(0  “  f'dampei'l  ^3  “  ^3  ^2 

111313 (t)  =  -*13(13  - 12)  -  C3  (^3  -  ij)  h  =  h 

=  2l  =  Xi-X 


1112X2(0  =  -  ^damper2  "I"  ^sprm^2  "I"  ^damper2 

1112X2(0  =  k3(X3  -  X2)  +  C3(X3  -  X2)  -*12 (X2  -  xj  -  “  Xi) 

1112(^2  +  Xi)  =  *13^3  +  C3Z3  -  k222  -  C2Z2 
m2Z2(t)  =  -m2Xi  +  ftgZj  4  C3Z3  -  ft2Z2  -  C2Z2 


3^lXi(0  -  “  ^dampers  "T  ^jpriJtpZ  "T  ^damper2 

ni^Xi  (0  =  ^2  (x2  -  Xj)  +  C2  (x2  -  x^)  -k|  (x^  -  x)  -  (x^  -  i) 
llll(%  +  x}  =  fe2Z2  +  C2Z2  - 


miZi(t]  =  -mixikiij  +  CjZi  -  k^Zi  -  qz^ 
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